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Abstract. We discuss the problem of a N two-level systems interacting with a single radiation mode in 
the strong-coupling regime. The thermodynamic properties of Dicke model are analyzed developing a 
perturbative expansion of the partition function in the high-temperature limit and we use this method to 
investigate the connections between the Dicke and the collective one-dimensional Ising model. 

PACS. 42. 50. Fx Cooperative phenomena in quantum optical systems - 05.70.Jk Critical point phenomena 
- 73.43.Nq Quantum phase transition 



1 Introduction 

The thermodynamic properties of the Dicke model [1] for 
N two-level system interacting with a single mode of the 
radiation field have been studied extensively by different 
authors and with different methods [2,3,4,5,6,7,8,9]. In- 
terest in this model has been renewed by a number of 
theoretical works, where it is discussed in connection with 
quantum chaos [10,11] and entanglement [12,13,14] and 
for various physical systems as photonic band gap mate- 
rials [15,16] and Josephson junctions [17]. 

The Dicke model exhibits a second-order phase tran- 
sition and, after the first derivation due to Hepp and Licb 
[2] , a simple computational method was provided by Wang 
and Hioe [3] based on the use of Glauber's coherent states 
[18] for the radiation field and on the assumption that, in 
the thermodynamic limit (N — > oo, V — *■ oo but p ~ N/V 
finite), the field operators can be treated as c- number func- 
tions. The Wang and Hioe method has been recently ap- 
plied by Lee and Johnson [16,19] to derive the thermo- 
dynamic properties for an extended version of the Dicke 
Hamiltonian incorporating spin-spin and spin-boson in- 
teractions. Perturbative methods have recently been pro- 
posed to perform thermodynamical calculations without 
to have recourse to the Wang and Hioe computational 
method. In Rcf. [20] it has been developed a perturbative 
expansion of partition function and a simple analytic so- 
lution is found for high coupling constant. The existence 
of this soluble model allows to show that the interaction 
of independent atomic spins with a resonant photon mode 
can be interpreted as an effective spin-spin interaction of 
long range nature. The similarity between the Dicke Model 
and the collective XY-model has been found in Ref. [14], 

Correspondence to: liberti@fis.unical.it 



using another perturbative approach in order to derive a 
temperature-dependent effective atomic hamiltonian. The 
prediction of a phase transition has been criticized as 
obtained for nonrealistic Hamiltonian, the simple Dicke 
Hamiltonian in which both the counter rotating terms and 
the diamagnetic term are truncated, that violated gauge 
invariance. As a matter of fact, the critical properties do 
not changes qualitatively when the counter rotating terms 
are taken into account [5] but, when the diamagnetic term 
is included, the Super-radiant Phase Transition (SPT) is 
forbidden to occur [21]. The same conclusion has been ob- 
tained in Ref. [20] for a model Hamiltonian incorporating 
the effects of the diamagnetic term into a new frequency 
of the photon mode [22]. 



In this paper we present a formal analysis which gen- 
eralizes the results of Ref. [20] beyond the resonant con- 
dition and gives an alternative derivation of some already 
known and recent results on the subject. In the next sec- 
tion we analyze the procedure that permit to develop a 
perturbative expansion of the partition function as gen- 
erally found in literature. The Dicke Model is introduced 
in Sec. 3 and a simple expression for the partition func- 
tion is analytically derived in Sec. 4 and 5 for high value 
of the coupling constant. This general result, as shown in 
App. A, is independent of the choice for the basis for ex- 
pressing the state of radiation field. Some special cases of 
interest are explored and the salient features of the ther- 
modynamic phase transition are discussed in Sec. 6 and 
in Sec. 7, where we also discuss the analogies between 
the Dicke model and the collective one-dimensional Ising 
model. In Sec. 8 we draw our final conclusions. 
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2 Perturbation expansion of partition function 3 The Model 



The starting point for deriving the thermodynamic prop- 
erties of a quantum model defined by an Hamiltonian H 
is the partition function: 



Z(N,T) = Tr{e- fiH } 



(1) 



It is generally difficult to diagonalize e~P H and a variety 
of systematic approximations has been proposed [23,24, 
27]. A convenient procedure is to introduce a separation 



H = H + Hi 



(2) 



such that the exponential operator of (1) can be disen- 
tangled into a product of an infinite series of exponential 
operators as 



-P(H +Hj) _ 



II' 



i=2 



(3) 



where Cj is a homogeneous polynomial of degree n in Hq 
and Hj. All the d s contain the commutator [Hi, Ho] 
and using the method given by Wilcox [26] they can be 
determined as: 

C 2 = \[Hj,Ho] (4) 
C 3 = —[Ho, [Hi, Ho]] - \[Hi, [Hi,H ]] (5) 

with increasing complexity for higher i. The basic step in 
the construction of approximants to the Eq. (3) is to find a 
product of exponential operators which is correct up to a 
certain power of (3. The disentangled and undiscntangled 
form of Eq. (3) are expanded in terms of /3 and operator 
coefficients of equal power of j3 are compared. The result 
is that 



0{H o +H!) _ ^pHo^DH, 



+ 0{(3 2 



(6) 



A more accurate approximation is obtained by introduc- 
ing the following symmetrized approximation of the Her- 
mitian operator e~@ H : 



-0{H o +H,) 



e -f]H a /2 e -f3H Ie -l3H /2 



+ 0(/3 3 ), (7) 



Obviously, the partition function obtained with using (7) 
and (6) is identical due to the cyclic permutation property 
of the trace. The error induced by the approximation 

Tr { e -/3(ffo+^) j = Tr{e-> 3H «e-> 3H '} (8) 
can be estimated using the Hermitian approximation 

e -/3(Ho+Hj) 

= e -f3H /2 e -pH I /2 e -p 3 C 3 /4 e -f3H I /2 e -f3H„/2 + 0(/3 5 )(9) 

and one obtains [24] 



| Tr | e -/3(H +H 7 ) _ e - / 3ffo/2 e -/3ff Je - / 3ff /2}| 



<£ |Tr{C 3 }| |Tr{e 



-t3H ap -PH, 



}| 



(10) 



We use this method to study the thermodynamic proper- 
ties of the Dicke model Hamiltonian (h = c = 1) 



N r 



H 



(11) 



Here, u> is the frequency of a single mode of radiation, e is 
the energy difference between the two levels of N identi- 
cal spin-1 systems, of, o~l and are respectively the z 
component, the raising and the lowering operators of the 
Pauli matrices used to describe the ith spin, a and are 
the annihilation and creation operators for photons. For 
N two-level atoms the coupling constant is 



A = ed 



2irp 



(12) 



where d is the projection of the transition dipole moment 
on the polarization vector of the field mode and p the 
density of the atoms. 

We make the following separation: 



Ho 
Hi 



where 



JS Z + -£=(at + a)(S+ + S~) 

2 y/N 



N 



S (*.±)=^ a (*. ± ) 



(13) 
(14) 

(15) 



i=l 



are the collective atomic operators. By applying the harmonic- 
oscillator commutation relations of the field mode opera- 
tors a T , a: 



[a, a t ] = l, [a^a, a*] = a) , [a) a, a] = — 
and that for the atomic operators: 

[S+,S-]=S Z , [S Z ,S ± ]=±2S ± 
it is easy to show that 



(16) 



(17) 



C 2 = —^(a - a.i)(S- +S+) 

[oj(a + a* ) (S+ + S-) + 2e(a - a f ) (S~ - S+) 

(18) 



C 3 = 



2VN 
4A 



+ —= (s-+s+) 



In order to evaluate the quality of the approximation we 
calculate the trace of C3 that requires summation over 
both the atomic and the field variables. A sum over atomic 
variables is well-suited to calculating the trace, yielding 

J2 ■■■ J2 (Si...S n \C 3 \S 1 ...S n ) = U 2 uj (19) 

Si=±l Sjv=±1 
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The expression (8) is correct up to the third order in (3 (or 
rather for {3 3 \ 2 lj < 1), i.e. is an appropriate description 
for the high-temperature limit. 

Although it is possible to derive higher-order approxi- 
mations of partition function in a systematic manner, the 
increasing complexity of Cj for higher i now complicates 
the expressions in such way that Z(N, T) is very difficult 
to evaluate analytically and requires numerical calcula- 
tions [27] . It is important to note that also a good choice of 
decomposing the Hamiltonian may affect the complexity 
of the calculation and the error induced by the approxima- 
tion (8). A different separation was studied in Ref. [20], in 
the simplest case of exact resonance between atom energy 
levels and frequency of radiation (u> = e) . The approxima- 
tions that one obtains from this choice is derived splitting 
off from (14) the eS z /2 term too and is, for this reason, 
less accurate. 



4 Partition function: Atomic variables 

In the form of Eq. (8), the partition function may be per- 
formed analytically. Writing out explicitly the trace over 
the atomic variables, the partition function is given by 

Z(N,T) = Tr F I e -0<"° ta ■■■ 

[ Si=±l S N =±1 

x (S u . . . , S N \e~ p ^^ hj \S u ..., S N ) | (20) 



where 



Noting that this operator has the property 

[huhj] =0, (i^j) 
from which it follows that 



N 



(21) 



(22) 



(23) 



we can reduce the partition function to the simpler form 
Z(N,T) = Tr F j e - /Wa 



J2 (S\e-^ Z+ ^ +a) ^\S) 



,s=±i 



N ' 



(24) 



where a x = <j + +<t . Expanding the exponential operators 
in a power series, we obtain 



0- 



2 k 



k=0 



(2k)\ 



4 + TV 



(a* + af 



2k + 



(25) 



where we have used the following Pauli matrices properties 
o* = o* = I , a*o* + a x a z = (26) 
Therefore, the sum of Eq.(24) is given by 

Es= ± i(SK^ + ^ (at+a) ^ is) 

_ 9 ^°° £.[e ! j.AV„t i 
- Z l^k=0 JW)\ [T + IV V + a > 

5 Partition function: Field variables 



(27) 



Using the Fock-state |n) for the photon field, the partition 
function is given by 



Z(N, T) = 2 N J2 ex P (~/3wn) 



71=0 



HE 



>2k 



(2k)\ 



J+N {a+a} 



N 



\n). (28) 



Since the operator (a) + a) commutes with itself, the 
power series appearing into the above equation can be 
written 



00 a2k r,2 \2 



E 

K 



(2k)\ 



N 



nE 

V.t=l fe,=0 



/3 2 



(2h)\ 



where K = k\ + ■ ■ ■ + fcjv- The matrix elements of photon 
operators are 

f]2q 

(n\(a^ +a) 2q \n) = —(n\e^ a + a ^\n} 



ri=0 



(Pi 
df] 2q 



e 2 L n (-r] 2 



(30) 



■q=0 



where L n {x) is the nth Laguerre polynomial. At this point 
we are able to write down the partition function (28) as 



K\ 



e\2(if-9) 



X \ 2q d 2 « , 



e ^J2e-^L n (- V 2 ) 



n=0 



(31) 



r/=0 



The sum over n is well known [28] to be given by 
f>-^% l (-, 2 ) = T ^exp^_l 

n=0 



> l_ e -/Ja, 

Substituting this into Eq. (31) one obtains 



0" _ I 



(32) 



Z{N,T) 



II E rofc.v ) E 



1 - e-P' 



Kl 



/e\2(if-!) / A 

V2J ITf 



2<I 



d 2 1 



2L- cot hH£ 



(33) 

r;=0 
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At this stage we use the following result 



drJ 2q 



ri=0 



(2q- l)!!coth 9 ^, q>0 (34) 



and the integral representation for the double factorial [28] 



(2<Z-1)M 



So, one has 



Z{N,T) = 



1 



2i+ 



— f 



dze 



.29 



(35) 



1 - e~P" Vi^f 



dze 4 



/ N oo 

iie 

\i=l fe,=0 







2k, 



(2h)l 



E 

q=0 



ql(K-q) 



I (I) 



2(K-q) ( X 



2q 



/2NJ 



coth 9 



fa 



1 - e-f^ V4^f 



dze 4 




A 2 z 2 fa 
■^ o (2fc)! L4 ' 2N C ° th 2 

which can be written in the final form 

1 If 00 

Z{N,T) = s = / dze- 

1 - e-^ v^F J-oo 



N 



(36) 



x < 2 cosh 




X 2 z 2 , //3w 
— — coth — — 
2AT V 2 




where z is an order parameter. This result is independent 
from the choice of different states as a basis for expressing 
the state of radiation field and the same result can be 
derived using the coherent states \a) for the photon field 
(see App. A). 



6 Phase transition: classical limit 

The integral (37) may be evaluated in the limit N — > oo 
by the steepest descent method. Writing the integral in 
term of a new variable x = z/VN we search the value x 
for which 



fix) 



+ In < 2 cosh 



„ e 2 A 2 x 2 , (fa 



is minimized. The minimum condition implies 



(38) 



fa 2 tanh 



tanh 



. e 2 X 2 x 2 , (fa 



fa" 




^cothf^ 
2 V 2 



(39) 



The existence of a nonzero solution of the above equa- 
tion means a phase transition. Through Eq. (39) we can 
compute a critical temperature and we find 



Pc = 



e tanh(/3 c w/2) 
2A 2 tanh(/3 c e/2) 



(40) 



The quantity (S z ) is of physical interest. We calculate it 
as follows: 

(S,) = -~hiZ(N,T) (41) 



From Eq. (39) we infer that 

- tanh (f ) , 13 < P c ; 
-2^tanh(^),/3>/3 c . 



N 



(42) 



A special case of interest is the limit reached when fa <C 
1. In this limit the partition function, which is given by 
Eq. (37), becomes 



Z{N,T) 



Att fa 



2 cosh (3 



dze 




N 



Let 



one has 



AN fa 



(43) 



(44) 



Z(N,T) ~ 



fair 



f 

J — < 



dye 



-Nf3ujy 2 



2 cosh | (3\j— +4X 2 y 2 



N 



(45) 



that correspond to the partition function obtained with 
the Wang and Hioe computational method. In this limit, 
the square of the order parameter y represent the average 
number of photons. Just as for the Wang and Hioe result 
[5] , the critical temperature is obtained from the equation 



tanh 



8 c e 
~2 



e lu 
4A 2 



(46) 



and the model Hamiltonian (11) undergoes a phase tran- 
sition at the critical value of the coupling constant A c = 
y/euj/2. For a coupling A > A c , above the critical tem- 
perature, the system is in the " normal phase" , whereas 
for P > (3 C , the equation (39) has a solution i / and 
the system is in the so-called "super-radiant phase". In 
this region, the expectation value of the collective angular 
momentum operator (42) is 



TV 



tanh ) , < p c ; 

0>0c- 



(47) 
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7 Phase transition: Ising limit 

The equivalence between the Dicke Model and the one- 
dimensional Ising Model for a system of mutually inter- 
acting spins 1/2 embedded in a transverse magnetic field 
[29] has been studied with different techniques [13,14,20]. 
This similarity emerges in our approximate scheme when 
/3e <C 1, where Eq. (40) becomes 



tanh 



If) = ^ 



(48) 



This result may be easily obtained observing that, in this 
limit, our model is equivalent to a high-temperature ex- 
pansion of the temperature-dependent effective Hamilto- 
nian given by 



H{p)=wa!a+-S x 



px 2 

2N 



coth 



0w 



(49) 



This effective Hamiltonian is the lowest-order effective 
Hamiltonian that can be obtained using the method of 
Ref. [14] (see App. B). As an alternative to the expression 
of Eq. (37), the partition function can be written as 

1 If 00 * 2 
Z(N,T) = 7 . = / dze~~ 

IA /coth(%i)i 



x Tr A { e 



(50) 



i.e., in the high temperature limit, 

1 1 f°° z 2 

Z(N,T) ~ 7 . = / dze-~ 

v ' i - e-v« 7_oo 

xTr A | e -^e^^ h (^)^| 



1-e- 



-Tr A { e-^ SZ e^ coth ^) Sl 



.51) 



As in the most general case (37), we can write the integral 
in term of x = z/y/N and utilize (in the limit N — > oo) 
the steepest descent method in order to search the value 
x for which 



/(*) = -Y + ln { 2cosh (y 



x cosh 



(3\z 



coth I — ^ 
2AM/ V 2 



(52) 



is minimized. The minimum condition implies 

[3\x 



x 
2 



' coth 



(3u 
~2~ 



tanh 



V2 



/ coth 



(3uj 
~2 



(53) 



The existence of a nonzero solution of the above equation 
means a phase transition. Through Eq. (53) we can com- 
pute a critical temperature that we find to be identical 



to the result (48). The order parameter near the critical 
temperature is 



2V3 



10- & 



0c 



(54) 



In the limit (3 c uj -C 1, we get (3 C ~ w/2A 2 and the interac- 
tion of independent spins with a photon mode induce an 
effective spin-spin interaction of long range nature that 
may be described by the Hamiltonian 



H = 



2l. 

Nlu' 



(55) 



that is in agreement with the results of Ref. [10,13]. 

Before leaving this section, we note that, for the special 
case of exact resonance u = e, Eq. (40) reduces to 



Pc ~ 2A 2 



(56) 



and this lead to the absence of a critical value of the cou- 
pling constant, i.e. to the result that the phase transition 
could occur even for A < e/2. However, the approxima- 
tion that we have proposed is valid for 0^.\ 2 uj < 1 and 
this condition permit us to derive accurate results only 
for A > e/2, i.e. for coupling sufficiently higher that A c . 



8 Conclusion 

The thermodynamic properties of a system governed by 
the Dicke Hamiltonian have been treated in the frame- 
work of an approximate model, through a perturbative 
expansion of partition function obtained by decomposing 
the Hamiltonian into two non-commuting hcrmitian oper- 
ators. This technique represent a practical and convenient 
method to determinate the behavior of the Dicke model 
in the strong-coupling regime. Wc have obtained a simple 
analytic expression for the partition function and the crit- 
ical temperature is easily determined. This technique was 
then used to explore some limiting cases and old and re- 
cent results are derived in an elementary and unified way. 
Our result extend those of Ref. [20] beyond the resonant 
interaction condition and support the recent arguments 
on the subject [13,14] concerning the similarity between 
the Dicke and collective one-dimensional Ising model. 



A Coherent states 

In terms of coherent states \a) [18], the partition function 
is given by 



Z(N,T) = 2 N j *£(a\e 

( 00 p2k 



.to ( 2fc ) ! 



N 



N 



H57) 
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where a\a) = a\a) and 



/ 



era 



a)(a\ = 1. 



(58) 



Using the over-completeness of the field coherent states 
(58) and the result 

(a|e-' Wa |7} = (a| 7 )e- Q *-'[ 1 - cx P(-' 3 ^] (59) 
Eq. (57) takes the form 

**7[1 — exp(— 



J 7T J IT 



(71 E 



2A- 



.fe=0 



(2k)\ 



c 2 \ 2 

7 + ^ flt + °> 2 



AT 



K)50) 



The partition function may be derived from Eq. (29) by 
computing the matrix element 



( 7 |(a t +a) 29 |a) 



d 2 i 
drfi 



( 7 |e"( ot +°>|a) 



r/=0 



( 7 | a ) ^ e 4 e vW+°) 



and by using 
One obtains 



dr) 2 " 

(a| 7 )| 2 = e -l Q -Tl 



(61) 
(62) 



)-'/,•, 



K 

E 

9=0 



A"! 



«!(Jf-g) 



— I ^I2 e -|a| 2 -|7l 2 +7*«+«*7^' 3aj +';(7*+«) ( 63 ) 

Recalling that the integration measure is defined to be 
given by 



e\W-o) ( X \ 2q d 2q nl 
e 2 



VN J dr] 2 i 



d 2 a dada* d(Rca) d(Ima) 

7T 2l 7T 



(64) 



one finds the partition function expression of Eq.(33), i.e. 
the same result obtained with the use of Fock states as a 
basis for the photon field. Finally, we want to underline 
that our results are derived without the c-number substi- 
tution for the field variables a — > a, the Wang and Hioc 
computational method, that permits to deal with non- 
interacting atoms subjected to an external magnetic field 
described by the amplitude a. 



B Effective Hamiltonian 

In this Appendix we will discuss the approach to the prob- 
lem of constructing effective atomic Hamiltonian for the 



Dicke Model. The Zassenhaus formula (3) can be used to 
obtain the following result: 



-0Hi 



i=2 



= E 



(-13)' 



i=0 



where 



fe=2 



(65) 



(66) 



The partition function of the whole system can be written 
as [27]: 



Z(N,T) = Tt a 



Tr F e 



-0H O 



E 

i=0 



(-/?)' 



= TrF (e-^°)TV A (e^ H '") 



(67) 



where 



H 



eff 



\i=0 



f*h q 



E- 



(Pi)F 



(68) 



where the thermal averaging is carried out with respect to 
the field variables, i.e. (0) F = Tr F ( e -P Ha O) /Tr F (e^ ). 
Eq. (68) can be rewritten as 



H 



j f = \ - P 9 1 



E 

3=1 



Q q 



(69) 



where the lowest-order terms arc 



Qi = (Pi)f = (H i ) f = -S z , 
Q2 = (Pi)l - (P2)F = (Hj) 2 F - (Hf) F 
= 4coth(^W 

Therefore the lowest-order effective Hamiltonian is 



2N 



H B /f = ls z -^cothl^.)Sl 



(3oj 
2 



(70) 



(71) 



in agreement with the result of Ref. [14] . 



C Inclusion of diamagnetic term 

In this Appendix we briefly discuss the inclusion of dia- 
magnetic effects in our model. Following the discussion 
of Ref. [22] and [20], the diamagnetic term may be in- 
corporated in the Dicke Hamiltonian of Eq. (11), without 
adding significant complication of the problem, by using 
a new field- mode of frequency fi instead of frequency ut, 
which is given by 



tt = sju(w + 4k) 



(72) 
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model Hamiltonian becomes 



N 



where tok — e 2 np/m. In terms of this new frequency, the 25. G. Polatsek and K.W. Becker, Phys. Rev. B 55, 16096 

(1997) 

26. R.M. Wilcox, J. Math. Phys. 8, 962 (1967) 

27. P. de Vries, A. Lagendijk and H. De Raedt, Phys. Rev. B 
(73) 35, 5425 (1987) 

28. G. Arfkcn, Mathematical Methods for Physicists (Aca- 
demic Press Inc., San Diego, 1985) 

29. G. Parisi, Statistical Field Theory, (Addison- Wesley P. C, 
Redwood, 1988) 



H = Qa)a + ^ 



|^ + --^(« t + «)K + + ^ ) 



where 



2itp 



(74) 



The Hamiltonian (73) is formally identical to the Hamil- 
tonian (11) and can be used to derive thermodynamic re- 
sults. The phase transition that is obtained in the limit of 
high coupling constant (A 3> A c = ^/U}/2) cannot occur 
due to sum-rule arguments [21] that requires A < 
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